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PREFACE 


Part  of  the  research  program  of  The  RAND  Corporation 
consists  of  basic  supporting  studies  in  mathematical  physics. 
One  area  of  this  is  concerned  with  the  study  of  turbulence. 

In  this  field,  the  analysis  of  "model  equations"  has  proved 

to  be  quite  useful. 

In  a  series  of  Memoranda,  of  which  this  is  the 
second,*  computational  techniques  which  approximate  the 

solution  to  the  model  equation  of  Burgers  are  investigated. 

I 

. 


*See  Ref.  2. 


The  nonlinearity  of  the  equations  of  turbulence 
forces  the  mathematician  to  study  approximation  techniques, 
such  as  "model  equations"  which  exhibit  many  of  the 
characteristics  of  the  realistic  equations. 

In  this  series,  the  authors  investigate  computational 
methods  £o  the  solution  of  the  model  equation  of  Burgers. 
In  the  present  Memorandum  the  application  of  a  difference 
algorithm  is  studied.  Numerical  results  obtained  from  a 
FORTRAN  program  are  given. 


ANOTHER  COMPUTATIONAL  APPROACH  TO  A  MATHEMATICAL 
MODEL  OF  TURBULENCE 


1.  INTRODUCTION 

In  the  study  of  turbulence  it  is  often  quite  useful 
to  analyze  '’model  equations"  which,  hopefully,  exhibit 
many  of  the  characteristics  of  the  more  realistic  equations. 
One  of  the  best  known  of  these  model  equations  is  that  of 
Burgers  which  is  given  by 


\  +  uux 


cu 


XX 

u(x,0)  -  g(x)  . 


(1.1) 


This  equation  is  discussed  by  Hopf  [13,  and  in  the 
first  paper  of  this  series  [2]  a  numerical  technique  is 
presented  in  which  this  equation  is  converted  into  an  in¬ 
finite  system  of  ordinary  differential  equations. 

In  this  second  paper,  another  numerical  technique  is 
investigated.  Here,  an  approximating  algorithm  is  used 
which  estimates  the  solution  of  (1.1)  to  an  accuracy  of  at 
least  0(A2),  where  A  is  the  integration  stepsize.  Similar 
methods  were  studied  in  previous  papers  [3,4],  and  shown 
to  be  quite  useful. 


Consider  the  equation  (1.1)  over  the  region  0  <  x  <  1, 
t  >  0,  and  suppose  that  g(x)  Is  periodic,  with  period  tt. 
Let  the  approximating  algorithm  be  given  by 


u(x,t+A)  -  Xu(x-au(x,t)4,t) 


+  [u(x+b » t)  +  u(x-b*\/A,t)] 


(2.1) 


where  A  Is  the  Integration  stepsize,  and  X,  a,  b  are  con¬ 
stants  which  will  be  determined.  To  show  that  (2.1) 
approximates  (1.1)  to  an  error  of  0(A2),  expand  both  sides 
of  (2.1)  In  a  Taylor  series  to  the  A2  term,  obtaining  the 
equation, 


ut  -  -Xauux  +  (1-X)  j-  u^ 


(2.2) 


For  (2.2)  to  approximate  (1.1),  the  following  relations 


must  hold: 


(2.3) 


(2.4) 


varied.  As  can  be  seen  in  the  following  table,  one  obtains 


essentially  the  same  results  for  M  ■  R  ■  10,  as  compared 
to  M  -  R  -  15,  in  less  than  half  the  time. 


nr 

t 

u(x 

ft 1 _ 

■ 

fll 

M-R-10 

M-R-15 

m 

-.355 

-.354 

19 

189 

-.319 

-.326 

.6 

3.00 

-.116  ! 

-.116 

.9 

4.00 

-.022 

-.022 

.2 

■SI 

-.132 

-.133 

Time 

40  sec 

1  min 

Note  that  in  this  example  and  in  the  following  examples 
the  largest  differences  occur  for  small  values  of  x. 


As  shown  in  Fig.  1,  varying  X  displays  large  differences 
in  the  values  of  u(x,t)  for  small  x  and  for  small  t.  Note, 

however,  that  the  maximum  of  each  curve  occurs  at  about 

. 

the  same  time.  For  small  t  and  large  x,  the  variation  in 
u(x,t)  is  much  smaller  (see  Fig.  2);  and  for  large  values 
of  t,  the  variation  is  quite  small  (see  the  table  below). 
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fl 
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.1 
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-.115 

.017 
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c)  u  +  uu  -  eu  where 

A  AA 

u(x,0)  ■  -sin  ttx ,  0  <  x  <  1 


c  -  .01 


V 


.  . 


The  parameter  which  has  the  greatest  effect  on  the 
results  is  the  integration  stepsize  A.  In  the  following 


table  the  differences  can  be  seen. 


4 .  HIGHER-ORDER  APPROXIMATION 

The  general  form  for  the  approximating  algorithm  is 
given  by 

N 

u(x,  t+A)  -  Xu  (x-au(x,t)A,t)  +  7  £u(x+b^’\/A,t) 
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.  ■< 


'  v' 
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in  which  the  X,  R,  «1#  and  a  arc  input*.  It  wa* 
shown  in  [5]  that  higher- ordered  approximation* ,  a*  in 
(4.1),  can  give  more-accurate  results,  provided  the  poly 
nooial  approximation  (2.6)  is  sufficiently  accurate. 
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